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Electrostatic force microscopy at cryogenic temperatures 
was used to probe the electrostatic interaction between a con- 
ductive atomic force microscopy tip and electronic charges 
trapped in an InAs quantum dot. Measurement of the self- 
oscillation frequency shift of the cantilever as a function of tip- 
sample bias voltage revealed discrete jumps at certain volt- 
ages when the tip was positioned above a quantum dot. These 
jumps are attributed to single-electron filling of the electronic 
states of the quantum dot. The estimated resonant energies 
of two s-states and four p-states are compared with the ex- 
perimental data obtained for an ensemble of quantum dots by 
the capacitance measurement technique. 

PACS numbers: 39, 68.37.Ps, 82.37.Gk 



I. INTRODUCTION 

Self assembled InAs quantum dots (QDs) grown on a 
GaAs substrate have attracted much attention recently 
due to their high crystal quality and numerous potential 
applications. Optical and electrical spectroscopy tech- 
niques to measure the discrete energy levels of InAs QDs 
have been developed Optical spectroscopy allows 

high energy resolution, and individual QDs can be stud- 
ied by spatially isolating them. This technique probes 
only the allowed transition of electron- hole pair. Capaci- 
tance spectroscopy, on the other hand, can probe electron 
and hole resonant energies seperately, but has so far been 
performed only on structures containing large numbers of 
quantum dots, due to the lack of spatial resolution and 
sensitivity. Recently, scanning capacitance microscopy in 
contact mode has been applied to capacitance imaging of 
QDs without evidence for single charge effects || . 

Application of atomic force microscopy ( AFM) to mea- 
surement of electrostatic force gradients due to few elec- 
tronic charges trapped below the surface has been demon- 
strated . In these experiments, only indirect esti- 
mations of the electron energy levels were given based 
on tunneling rate from the localized states into the 
substrate. Measurement of frequency shift of a self- 
oscillating cantilever as a function of bias voltage has also 
been used to look at sub-surface charge densities without 
indication of single charge effects || . 

In this letter, we report, a new measurement technique 
to detect trapped charge in individual InAs QDs. A non- 
contact measurement scheme is employed, where a con- 
ductive AFM tip is positioned above a QD to probe the 
gradient of the electrostatic force between the tip and 



charges trapped in the QD as a function of tip-sample 
potential difference. 



II. EXPERIMENT 

The sample used in the experiment was grown on an 
n-type GaAs (100) substrate by MBE. A 0.5 ^m thick 
GaAs buffer layer with 10 18 cm -3 Si doping was first 
grown, followed by the growth of a semi-insulating GaAs 
tunnel barrier of 15 nm thickness. After growth of a 
monolayer InAs wetting layer, InAs QDs with average 
base diameter of 20 nm were grown on top of the tun- 
nel barrier. The density of the QDs was observed to be 
about 10 10 cm~ 2 at the center of the wafer and decreas- 
ing towards the edges. The homemade AFM used in this 
experiment used a fiber-interferometric deflection sensor 
for cantilever readout, with a wavelength of 1310 nm and 
an incident light power of about 5 /iW. 
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FIG. 1. (a) Contact mode image of InAs quantum dots. 
Locations where electrostatic force spectroscopy is performed 
are labeled A and B. (b) Schematic of the experimental setup. 
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The Pt/Ir coated AFM cantilever had a spring con- 
stant of 2.5 N/m and a fundamental resonance frequency 
of 68.8 KHz. The Q factor of the cantilever was about 
28,000 at 77K and increased to 59,000 at 4.2 K. Contact 
mode images of QDs were obtained prior to electrostatic 
force spectroscopy (EFS) on individual dots. Imaging 
and EFS were done at a temperature of 4.2 K. Since 
there was no capping layer on the QDs, base diameters 
and heights can be accurately determined from the con- 
tact mode images (Fig. la). 

After contact mode imaging, the tip was lifted 10 to 
30 nm above the surface. Self oscillation of the cantilever 
was sustained through an analog tracking oscillator cir- 
cuit ||. The circuit is a phase-locked loop (PLL) sys- 
tem that drives the cantilever at its resonant frequency 
with a fixed drive amplitude. Tracking of the frequency 
is accomplished by locking the local oscillator phase to 
the AFM cantilever oscillation phase by a proportional- 
integral feedback controller. The block diagram of the 
EFS setup is shown in Fig. lb. During the EFS measure- 
ments, the oscillation amplitude was fixed to be about 1.5 
nm, to provide adequate signal-to-noise ratio while main- 
taining a small ratio between the oscillation amplitude 
and the tip-sample seperation. In this regime, the lin- 
ear approximation used in the theoretical calculation of 
the frequency shifts is justified. The minimum detectable 
frequency shift is fundamentally limited by the thermo- 
mechanical fluctuation of the cantilever and is given by 
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where fo is the cantilever resonance frequency, k B T is 
the thermal energy, Q is the mechanical quality factor, k 
is the spring constant, A osc is the amplitude of cantilever 
oscillation and B is the measurement bandwidth. The 
detectable frequency shift was also limited by the noise 
sources in the secondary detectors such as optical inter- 
ferometer and electronic amplifier noise, and the overall 
sensitivity of our system was about 0.5 Hz/VHz. 

Since the tip radius is about 20 nm and the EFS is 
performed at a tip sample separation of about 15 nm on 
average, a parallel plate capacitor model can be used to 
analyze the experimental data. The force between the 
tip and sample in the presence of a trapped charge Q in 
the QD is given by (Refs. §§): 
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Here, A is the effective tip area, di is the tunnel barrier 
thickness, d 2 is the optional capping layer thickness, V ts 
is the tip-sample bias voltage, e r is the relative dielectric 
constant of a GaAs tunnel barrier, and z is the tip to 
sample-surface separation. A WKB calculation based on 
a bulk and QD band diagram suggests that the tunneling 



rate into the QD is fast enough that the QD is charged 
and discharged at a time scale much faster than other 
time scales involved in the experiment. Therefore, we as- 
sume that as the tip-sample voltage is swept, the dot and 
the n-type reservoir are in equilibrium with each other. 

Ins pec ting individual terms of the electrostatic force, 
Eq. (2.2), it is seen that only the first and second terms 
of the sum contains information about the localized state 
and the third term merely provides a background elec- 
trostatic force from the substrate. For experimental con- 
ditions of interest, it is also seen that the second term of 
the sum is two orders of magnitude larger than the first 
term, therefore the first term can be neglected. Inclu ding 
only the dominant term, we can approximate Eq. ( |2.2| ) 
for a number of trap states that are at distances hj above 
the substrate with individual charges qj as 
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where d is the overall dielectric thickness covering the 
n + substrate. To calculate the trapped charge q^ for a 
given state i as a function of tip-sample bias voltage, it is 
assumed that the quantum dot is in equilibrium with the 
substrate electron distribution. For state i, the charge in 
the state can be calculated by the Fermi-Dirac distribu- 
tion 



= cfl 



exp( 



E f 



k R T 



)]' 



(2.4) 



In the parallel plate approximation, ignoring band 
bending effects at the interfaces, energy level Ej can be 
expressed with reference to the bulk fermi level as the 
sum of a bias induced energy shift and an inherent en- 
ergy Ej . Assuming the bulk fermi level Ef to be the 
conduction band energy Ecb for the n + substrate, this 
bias dependent energy for state i is given by 



E i (V tB )=E i0 +V t 
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As can be seen in Eq. (|2.5D , Ej is the energy of state 
i with respect to the bulk conduction band Ecb under 
the zero bias condition. The voltage dividing ratio which 
multiplies V ts in Eq. (2.5) can be deduced from the thick- 
nesses of the barrier layers and tip-sample separation. In 
our experiments this ratio was typically 0.03 to 0.1. 

During the course of EFS, the resonance frequency of 
the cantilever is shifted in the presence of the external 
electrostatic force gradient. This shift Jf can be calcu- 
lated in the limit of small oscillation amplitude as 
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One can approximate this frequency shift 5i as 
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where the bracket denotes the average of the electro- 
static force gradient over the oscillation amplitude of the 
cantilever. 

In order to isolate the discrete charging signature, we 
look at the derivative of the frequency shift with respect 
to the tip sample voltage 



d 



■Si: 



fo 



2k A o 



dV t 



■[F e (z + A osc /2) -F e (z- A osc /2)} 



(2.8) 



To calculate the con tribution of electron state i t o the 
derivative in Eq. (2.8), inserting F c j(z) from Eq. (2.5) 



and neglecting small terms, we get 
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(2.9) 



In the zero temperature limit, Fermi-Dirac distribu- 
tion which determines the charging of the electron states 
approaches a step function. Therefore, in this limit, Eq. 
(2.9) predicts two delta functions for each state i below 
the tip, appearing at tip-sample voltages that depend on 
tip-sample separation, overall semi-insulating layer thick- 
ness, the oscillation amplitude of the cantilever as well as 
the spatial location of the state and its energy with re- 
spect to the bulk fermi level of the n + GaAs region. The 
peaks in d(6£)/ffVts occur when the condition 
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is satisfied. 

In the finite temperat ure limit, further information can 
be obtained from Eq. (2.9) as it estimates two peaks of 
opposite sign with magnitudes and widths determined 
by the height hj of level i above the n + substrate. The 
derivative <9qj(z, Vt s )/9V ts is 
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by 



The width of the individual peaks can be approximated 
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By analyzing EFS data for different tip-sample surface 
separations and cantilever oscillation amplitudes in a pa- 
rameter range where the parallel plate approximation is 



valid, using Eq. ( 2.1 0| ) and Eq. (2.12) one can solve for 



hj and Ej and extract information on the location and 
the energy of states. Repeating the procedure at multiple 
locations on the sample, a 3D mapping of the localized 
density of states can be obtained. 
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FIG. 2. Electrostatic force spectroscopy at (a) location A 
and (b) location B of Fig. la . The numerical derivative 
of frequency shift with respect to bias voltage, d(<5f)/dV is 
plotted versus bias voltage to clarify jumps in the cantilever 
resonance frequency. Insets of (a) and (b) show the actual 
frequency shift data on the quantum dot and on the wetting 
layer. 

The contact mode image shown in Fig. la shows two 
locations where the EFS was performed at a tip-sample 
surface separation of 13 nm. The obtained frequency 
shift data are plotted in Fig. 2. Calibration of the tilt 
of the sample prior to EFS measurements ensured that 
the tip-sample separation were identical for these two 
points. The drift of the scanner was also measured to be 
insignificant during the data acquisition period. Discrete 
jumps in the cantilever frequency were observed on top 
of a QD having a 40 nm base diameter and a height of 
17 nm. These jumps disappeared as the tip moved away 
from the QD in the plane of the surface, while the height 
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was kept constant. 




FIG. 3. The derivative of frequency shift with respect to 
bias voltage, d(5f)/dV calculated from the theory described 
in the text. Tip-sample separation and oscillation amplitude 
used in calculation are 14.2 nm and 1 nm respectively. Tip 
radius is assumed to be 20 nm. Electron energy levels are 
obtained from the data presented in Table I . 

Based on the calculations of QD electron energy levels, 
we expect to see two s-like states whose energy degen- 
eracy is lifted by an on-site Coulomb charging energy. 
In addition, four p-like states whose energy degenera- 
cies are lifted both by on-site Coulomb charging effect 
and strain-induced anisotropy are expected. According 
to calculations done for various sized QDs JlCj , the on- 
site Coulomb charging energies for s-like and p-like states 
differ and are estimated to be about 18 meV and 10 meV 
respectively, for 20 nm base diameter QDs in a semi- 
insulating GaAs matrix. 

TABLE I. Electron energy levels inferred from previous ca- 
pacitance-voltage (CV) spectroscopy for ensemble of 20 nm 
average base diameter capped QDs and EFS experimental 
data presented in this paper involving 40 nm base diameter 
uncapped QD. Electron energies are in units of meV and are 
shifted to match the ground state energies, E s i- 
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The charging energies decrease to 11 meV and 5 meV 
for the s-like and p-like states in 40 nm base diameter 
QDs. A first-principles calculation jll| that takes into 
account the strain anisotropy in the QDs predicts the 
four p-like states to form two pairs of states that are 
separated by about 25 meV. 

Applying the above mentioned model to the experi- 
mental data shown in Fig. 2, allows estimation of elec- 
tron energy levels for the first two states (s-like states) 
of 148 and 157 meV with respect to the InAs conduction 
band minimum. The four p-like states are also identified, 
and calculated to be 188, 194 , 216, and 222 meV above 
the InAs conduction band minimum. Since band bend- 
ing effects have not been considered in the calculation, 
and the work function difference of the tip and sample 
is known only approximately, these energies are correct 
only to within an overall offset. Using these QD energies 
in the electrostatic model described above, a theoretical 
calculation of the expected frequency shift versus bias 
voltage is plotted in Fig. 3. The data in Fig. 2a were 
measured on top of a QD having a 40 nm base diameter. 
The observed charging energies of 9 meV for the s-like 
states and 6 mev for the p-like states in EFS data agrees 
well with theory in Ref. 10. A comparison of electron 
energy levels for 20 nm base diameter InAs QDs mea- 
sured through conventional capacitance spectroscopy on 
ensembles of QDs (Ref. 4) with results of EFS measure- 
ment on a single 40 nm base QD is given in Table |. 

III. CONCLUSION 

In conclusion, we have demonstrated an electrostatic 
force microscopy technique for identifying the charge 
state of a single QD. Our results on uncapped dots of 
40 nm base diameter agree well with the theoretical pre- 
diction for the charging energies for the s-like and p-like 
states and and agree qualitatively with theoretical esti- 
mates of quantized energy levels for this size QDs. The 
EFS method presented has the advantage that it can 
be applied without ultra-high vacuum (UHV) conditions 
since the states probed can be below a capping layer. It 
can also be extended to the study of states at interfaces 
of properly designed MBE structures and of surfaces with 
adsorbates. 
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